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We derive analytical forms for non-flow contributions from cluster correlation to two-particle 
elliptic flow (i>2{2}) measure. We also derive an analytical form for jet-correlation flow-background 
with the same cluster approach. We argue that the elliptic flow vi parameter to be used in jet- 
correlation background is that from two-particle method excluding non-flow correlations unrelated 
to the reaction plane, but including cross-terms between cluster correlation and cluster flow. We 
verify our result with Monte Carlo simulations. We discuss how one may obtain the V2 parameter 
for jet-correlation background experimentally. 



PACS numbers: 

I. INTRODUCTION 

Jet-like angular correlation studies with high transverse momentum (jpt) trigger particles have provided valuable 
information on the properties of the medium created in relativistic heavy-ion collisions [H, 0. In such studies, 
correlation functions are formed in azimuthal angle difference between an associated particle and a high px trigger 
particle, which preferentially selects (di-)jet. One important aspect of these studies is the subtraction of combinatorial 
background which itself is non-uniform due to anisotropic particle distribution with respect to the reaction plane- 
both the trigger particle and the associated particles are correlated with the common reaction plane in an event. The 
critical part is to determine flow parameters, mainly elliptic flow (1)2), to be used in constructing background. 

There are many w 2 measurements [1, H[. They contain various degrees of non-flow contributions, such as those from 
resonance decays and jet correlations. Those non-flow effects should not be included in jet-correlation background. We 
shall refer to this jet-correlation background as flow-background. The anisotropic flow to be used for flow-background 
should be ideally that from two-particle method, t>2{2} 0,0, because jet-like correlation is analyzed by two-particle 
correlation method. Moreover, two-particle anisotropic flow contains fluctuations which should be included in jet- 
correlation flow-background 0, @ . 

Non-flow is due to azimuthal correlations unrelated to the reaction plane, such as resonances, (mini)jets, or generally, 
clusters. In this proceedings, we study non-flow contributions in two-particle i>a{2} in a cluster approach as in [7|; 
analytical form is derived for each non-flow component. We shall demonstrate that the flow to be used in jet- 
correlation background subtraction should be the two-particle t>2{2} excluding cluster correlations unrelated to the 
reaction plane, but including cross-terms between cluster correlation and cluster flow p|. We verify our analytical 
result with Monte Carlo simulations. We discuss how one may obtain the elliptic flow t>2 parameter for jet-correlation 
background experimentally. 

II. NON-FLOW EFFECT FROM CLUSTER CORRELATIONS 

Suppose an event is composed of particles from hydro-medium and clusters of various types (such as minijets and 
resonance decays). Particle pairs can be decomposed into four sources: 

• particle pairs from hydro-medium (B), 

particle pairs from same cluster (C), 

particle pairs between hydro- medium and clusters (X), and 
particle pairs between clusters (Y). 
The total sum of the cosines of pair opening angles is 



where 



^]cos2A0y=B+ C + 2X+ F ' W 

i^j k^cluster k^cluster (fci ) ^cluster 



B = ^2 cos2A4>ij, (2) 

(ijtj)ehydro 
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C = cos2A^, (3) 

X = Y cos2A ^> ( 4 ) 

i£k j^hydro 

Y = cos2A0 tf . (5) 

iefei jefe2 

Here i,j are particle indices, Afaj = (f>i — (f>j, and fc stands for a cluster. Below we derive analytical form for each 
source. 

A. Background flow correlation 

Hydro-background particle correlation is only from hydrodynamic anisotropic flow: 

B= Y cos2A(j)ij =V h y(cos2A<p h y) =VhyV%{2} hy (6) 

where Vhy = (Nh y (Nh y — 1)) is the number of background pairs. 

B. Particle correlation within cluster 

Particle correlation within cluster is given by 

r 2Tt 



C= V cos2A0jj = / V a {4>k)Pci((f>k)d(t>k x 



/ / o (A&,0 fc ) A & / f a {A<p j ,fa)dA<p j cos2(Acj> i -Acj> j ). (?) 
Jo Jo 

Here A(f>i,j = — 0fe is the azimuthal angles of particles in cluster k relative to cluster axis (f>k (which can be 
defined just for convenience); f a (A(f), <pk) is the correlation function of (associated) particles inside cluster k relative 
to the cluster axis <f>k, generally dependent of the cluster axis <j)k = <t>k — V'flP relative to the reaction plane, and 
| ' f a (A(j), 4>k)A(j) = 1; V a (4>k) is number of (associated) particle pairs in cluster k, generally dependent of the cluster 
axis; p c i(4>k) is the density function of cluster k relative to the reaction plane, which we will assume is given by elliptic 
flow of clusters, and J Q V p c i{4>k)d4>k = 1- Note, there can be many types of clusters (e.g. jet-correlation, resonance 
decays); the subscript l cV stands for one type of clusters and we have omitted summation over all types of clusters 
from the formulism; in this work we will discuss only one type of clusters at a time. 
In general, 

C = V a (cos 2A<f>ij) cl (8) 

where V a is average number of pairs per cluster and (cos 2A0 ij ) c ; is the average cosine of twice pair opening angle in 
the cluster. 

If particles inside cluster are independent of each other except all of them are correlated with the cluster axis, then 
we can factorize the correlation terms and obtain 

C = J o 'Pa{k)Pd{k)dk U / o (A0,0 fc )cos2A0dA0j +lj / a (A0,0 fe )sin2A0dA^J 

= f P„(^ fc )(<cos2A0)| fc + <siii2A0)| fc )/> c ,(^ fc )d0fc. (9) 

Here (cos2A</>) < ^ and (sin2A0)^ fc (A<p = 4> — <j>k) ar e averages within cluster k, and are generally dependent of the 
cluster axis (frk- 

If the cluster correlation function f a (A(p,<pk) is symmetric about Acf) = 0, then (sin2A0)^ = 0. Further, in 
the special case where particle correlation in clusters does not vary with cluster location cf>k, i.e., V a = const, and 
f a (A(j>,4>k) = const., then 

C = V a (cos2A<f>} 2 cl . (10) 
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C. Background-cluster correlation 

Correlation between cluster particles and hydro-medium particles is given by 

p2n p2ir 
X COS2A0jj = / N hy ((f> k=ly 2,...)phy(4>hy)d<j>hy / Pcl(4>k)d(f>k X 

i€k jehydro J J 

N a (4> k )fa(^i,4>k)dA^i [cos 2(0, - (j> hy )] . (11) 



Here <phy = 4>hy — ipRP, and Phy(4>hy) is the density function of hydro particles relative to the reaction plane (i.e., 
anisotropic hydro flow); N a {<fi k ) is number of (associated) particles in cluster and is generally dependent of the cluster 
axis 4>k- For generality we have taken the number of hydro-medium particles Nh y (4>k=i,2...) to depend on positions of 
all clusters. Such dependence can arise in real data analysis, such as jet-correlation analysis, from interplays between 
centrality cut and biases due to selection of specific clusters. Rewriting <pi — (f>h y = Acf>i + 4>k — 4>hy, we have 

p2rr 

X = I Nhy(4>k=l,2,...)phy(,(t>hy) COS 2 4>hyd(f>hy X 







2tt p2ir 

Pci(4>k)d4> k / JV a (&)/ o (A& J 0fc)o»2(A0 j +&)dA&. (12) 



io Jo 

/•2-K n 27T 



i*Z1Z I>Z7T 

Here we have used / d(f> k p c i(cf> k ) / dA(j> k f a (A(j> k , (f> k ) sm2(A(f) k + <f> k ) = because of symmetries f a (A(j> k ,<f> k ) 
J° Jo 

f a (—A<f> k ,—<f) k ) and p c i{4>k) = Pci{~4>k)- Note, due to elliptic flow of clusters, cluster particles acquire elliptic flow 



2 /-27T 

V2,a = (cos2(0 - ip R p)) = — / Pci(4>k)d(f)k N a (4> k )f a (A(t) i ,(i> k )cos2(A(j3 i + <f> k )dA(f) i . (13) 

iV a Jo JO 

Using the notation in Eq. (|T3|) . we have 

X = N hy N a v 2 ,hyV 2 ,a = N hy N a v 2 {2} hy v 2 {2} a . (14) 

Here the product of the v 2 s includes flow fluctuation, and equals to the product of two-particle U2's. This is because 
v 2 {2\ of hydro-particles and cluster particles contain only fluctuation; non-flow does not exist between hydro-particles, 
nor between particles from different clusters. (Note, two 'clusters' can originate from a common ancestor, such as jet 
fragmentation into two p mesons which in turn decay into two pairs of pions. In our formulism, such 'clusters' are 
considered to be parts of a single cluster rather than two p-decay clusters.) 

Again, in the special case where particle correlation in clusters does not vary with cluster location <p k (jV„ = const., 
N hy = const., and f a (A(f),$ k ) = f a (A<j>)), Eq. (JTSJ) becomes 

v 2 , a = I pci{4> k )d(f) k cos24> k I f a {A(j)i) cos 2A0 i dA<fo = v 2 , c z(cos2A<?(>) c/ , (15) 

and we have 



o Jo 



X = N hy N a v 2 {2} hy v 2 {2} c i(cos2A(f>} c i. (16) 

D. Particle correlation between clusters 

Correlation between particles from different clusters is given by 
Y = ^2 cos 2A ^i 

p2n n2ix 

= / Pd{4>kx)d^ kl I N a (4> kl )f a (A<f>i,4> kl )dA(f>iX 
Jo Jo 

/ p c S k2 )d4> k2 N a (4> k2 )f a (A^> j ,4> k2 )dA(l )j [cos2(^ i -(f> j )}, (17) 
Jo Jo 
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where k\ and k 2 stand for two clusters. Rewriting 0, — cj>j = A<pi + (p^ — A<pj — <f>k 2 , we obtain 

Y = N 2 a vl a = N 2 a v 2 {2} a , (18) 

where v 2 , a is given by Eq. (fT3|) . Again the cluster particle elliptic flow squared in Eq. (fT8|) contains flow fluctuation. 
In the special case where particle correlation in clusters does not vary with cluster location <pk, we have 

y = A 2 « 2 {2} ci (cos2A0) 2 ; . (19) 



E. Summary of non-flow effect from cluster correlations 

To summarize, let us now obtain the relationship between two-particle elliptic flow v 2 {2} that is affected by non-flow, 
and the real hydro-type two-particle elliptic flow v 2 {2}} ly . Assuming Poisson statistics, Eq. {T]) gives 

N 2 v 2 2 {2} = N 2 y v 2 2 {2} hy + N cl N 2 (cos 2A</> ij ) cl + 2N hy N cl N a v 2 {2} hy v 2 {2} a + 
N cl (N cl - 1)A 2 V 2 {2} Q 
= (N hy v 2 {2} hy + N cl N a v 2 {2} a ) 2 +N cl N 2 ((cos2A& 3 ) cl -vl{2} a ) (20) 

where A = N} ly + N c iN a , N c i is average number of clusters, and we have taken distributions of total multiplicity and 
number of particles per cluster to be Poisson, so that (A(A — 1)) = A 2 and (N a (N a — 1)) = A 2 . We have taken 
the number of cluster pairs to be A c z(A c ; — 1) (i.e., not Poisson) so that the total number of pairs adds up to A 2 . 
Rearranging, we have 



m = (^v 2 {2} hy + ^ V2 {2} a )\ ^gl «cos2A^) ci - , 2 {2} a ) . (21) 



N "^" inv 1 N " z ^ ,a J 1 N 2 
For many cluster types, Eq. (|2"Tj) is generalized to 

(^M2}hy + ^^v 2 {2} a ) + ^^((cos2A^) c; - V 2 {2} a ). (22) 

V cl / cl 

We shall focus on the special case where all clusters are of the same type and particle correlation in clusters does 
not vary with cluster axis relative to the reaction plane. Using Eq. (fT5|) . Eq. (f2"Tj) becomes 

«2{2} = ^« 2 {2}^ + ^^ W2 {2} ci (cos2A0) d X 
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NclN 2 
A 2 



((cos2A^) c/ - v 2 {2} c; (cos2A< ? (>) 2 / ) . (23) 



Eq. ([23]) can be rewritten into 



v 2 {2} = v 2 {2} hy + 2 N ^ V Ncl ^ a v 2 {2} hy {v 2 {2} cl {co S 2Acj>) cl - v 2 {2} hy ) + 
(^) 2 («2 2 {2} c /(cos2A0) 2 ; - v 2 {2} hy ) + 

^0- «coB2A^)d " « 2 2 {2} c /(cos2A0) 2 ; ) . (24) 

The second term on the r.h.s. is non-flow (beyond hydro-flow) due to correlation between hydro-particles and clus- 
ter particles in excess of that between two hydro-particles, and the third term is that due to correlation between 
particles from different clusters. These non-flow contributions, which are beyond hydro-flow, can be positive or 
negative, depending on the relative magnitudes of background particle flow and cluster flow diluted by particle 
spread inside cluster. The non-flow contributions are positive when « 2 {2} c ;(cos 2A<fi) c i > v 2 {2}h y and negative when 
v 2 {2} c i(cos2A(j)) ci < v 2 {2\hy- This can be easily understood because if v 2 {2} c i = v 2 {2}h y , then the angular smearing 
of particles inside each cluster, (cos2A0) c ;, makes the angular variation of cluster particles less than that of hydro- 
particles, resulting in a negative non-flow contribution. If the net effect of cluster anisotropy and particle distribution 
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inside clusters, f2{2} c ;(cos 2A0) C /, equals to hydro anisotropy, then cluster particles and hydro-particles have the same 
angular variation relative to the reaction plane, resulting in zero non-flow from cross-pairs between hydro-particles 
and cluster particles and between particles from different clusters. 

The second part of the last term of Eq. ([24]) r.h.s., ^j^-v 2 {2} c i(cos2A<j)} 2 l , arises from the fact that number of 
clusters is fixed in order to have Poisson statistics for hydro particle multiplicity, particle multiplicity in clusters, 
and total number of particles @- It can be safely neglected because generally v 2 {2} c i << 1. The first part of the 

last term of Eq. (|24|) r.h.s., JV ^ a (cos 2A<^j) e ;, is non-flow due to correlation between particles in the same cluster. 
This non-flow contribution can also be positive or negative. If particle emissions within clusters are independent, 
(cos 2 A0ij ■) d — (cos2A0)^ , then Eq. becomes 

v 2 2 {2} = v 2 {2} hy + 2 ^ Ncl ^ a v 2 {2} hy (^{2} d (cos2A</») d - v 2 {2} hy ) + 

N 
N 2 [ 



(t»|{2} c/ (cos2A0)2 ; - v 2 2 {2} hy ) + 
1 - v\{2} cl ) (cos2A^. (25) 



In this case the non-flow contribution due to particle correlation within clusters can only be positive. 

We note that the non-flow contributions from the second and third term of Eq. ([24| r.h.s. have the identical 
azimuthal shape relative to the reaction plane as that of hydro-flow, because they arise from the common correlation 
of clusters and hydro-particles to the reaction plane. As a result these non-flow contributions will unlikely be separated 
from medium hydro-flow in inclusive measurement of azimuthal correlation. To separate these two contributions, one 
needs to identify clusters and measure two-cluster azimuthal correlation. In fact, elliptic flow is often defined as the 
second harmonic of particle distribution relative to the reaction plane, v 2 {RP} = (cos2(</>— tppp)} . For events 
composed of hydro-particles and clusters, we have 

v 2 {RP} = ^ (cos 2(0 - ip R p))hy + NcL j^ a (cos 2(0 - lp RP ))ci 

= ^v 2 {RP} hy + ^j^v 2 {RP} cl (cos2A(b) cl . (26) 

This is analogous to the terms in the first pair of parentheses on Eq. (|23[) r.h.s. except the latter contains flow 
fluctuation. The elliptic flow definition by Eq. (I26[) contains cluster contribution through angular spread of particles 
in clusters, (cos2A0) c ;, and anisotropy of the clusters themselves, v 2 {RP} c i. This raises question to comparisons 
often made between elliptic flow measurements and hydro calculations which may include flow fluctuation but does 
not include cluster correlations. 

III. ELLIPTIC FLOW FOR JET-CORRELATION BACKGROUND 

In this section, we derive an analytical form for flow-background to jet-correlation in the cluster approach, as 
used in our non-flow study above and in [7j. We suppose a relativistic heavy-ion collision event is composed of 
hydrodynamic medium particles, jet-correlated particles, and particles correlated via clusters. Hydro-particles, high 
Pt trigger particles, and clusters are distributed relative to the reaction plane (tp) by 

dN N 

— = — [l + 2« 2 cos 2(<£ -</>)] (27) 

with the corresponding elliptic flow parameter v 2 and multiplicity N. Particle azimuthal distribution with respect to 
a trigger particle is 

N t dA(f> dA(f> ^ dA6 ^ dAcj) dAcj) [ ' 

k^jetdclus k£jet 

where Acf> — tfi — cf> t . In Eq. (|28|) . dN a j et / 'dAcf> is jet-correlation signal of interest. All other terms are backgrounds. If 
trigger particle multiplicity is Poisson and effects due to interplay between collision centrality selection (usually via 
multiplicity) and trigger bias are negligible, then the background event of a triggered (di-)jet should be identical to 
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any inclusive event, without requiring a high px trigger particle, but with all other event selection requirements as 
same as for triggered events [Io| . Thus, we can use inclusive events to obtain flow-background: 

J_dN_ = ( dN hy ^ dNa^k, V dNaA dN a . ]et 

N t dAcj) ~ a dA(f> L> dA6 ^ dAcj) + dAcj) 1 ' 

\ k^jetdclus kCziet / . 

\ ' me 

where a is a normalization factor, often determined by the assumption of ZYAM or ZYAl (zero jet-correlated yield 
at minimum or at Acj) = 1) (l5j . and is approximately unity. The background is 



dN bg _ dN hy \ - dN a , k x - dN a ,k _ dN hy v^^ dN a 
dA6 ~ dAcj) ^ dA6 ^ dAcj) ~ dAcj) ^ cl dA(j> 

k^jetdclus kdjet cl 



where we have eliminated subscript 'inc' to lighten notation. We have summed over all cluster types l cP including 
jet-correlation, where N c i is number of clusters of type l cl\ Different cluster types include jet and minijet correlations, 
resonance decays, etc. 

The hydro-background is simply 

^ = ^1 (l + 2v 2 ,m, hy cos 2A0) (31) 
dAtp Zir 

where V2,t is elliptic flow parameter of trigger particles and V2,hy is that of hydro-medium particles. 
The cluster particles background is given by 

dN f 2 * { 2lx ~ ~ f 2Tr 1 

d<p tPt (<f> t ) d<f> kPcl (<f> k ) dA(pif a (A(j)i, <p k ) x —6(A(f>i + cf> k — A<f) - cp t ) (32) 



dA(p J Q J J Q 2tt 

where cj> t = 4>t - ip, 4>k = 4>k - ip, Afa = fa - cj> k , and p t {fa) = ^ U + 2v 2 ,t cos 2cj>^j and p c i(<j>k) = 

^1 + 2v2. d cos 2<fikj are density profiles (i.e., ^-modulated distributions) of trigger particles and clusters rela- 
tive to the reaction plane, respectively. We have assumed that the cluster axis (or cluster parent) distribution is also 
anisotropic with respect to the reaction plane. In Eq. (|32|) . f a {Acj>i,cf> k ) = d J^ c j k is distribution of particles in cluster 
relative to cluster axis (cluster correlation function), which may depend on the cluster axis relative to the reaction 
plane cj) k Decomposing pt{<f>t), we obtain 

dN n 1 f 2v ~ ~ f 2 * 

'hPd(4>k) dAfaf a (dAfa,(f> k ) + 



dAcj) 2n Jo 



kpd(4>k) dAfaf a (Afa,rj) k )cos2(Acj) l + ci k - Acf>). (33) 



2vr 

Because of symmetry, f a (Afa, (j> k ) = f a (-Afa, -cj> k ) and p c i{cj) k ) = Pd(-<Pk), we have 



kPd{4>k) dAfaf a (Afa,4, k )sm2(Afa + fa) = 0. (34) 



Therefore 



dN n , 1 ~ 



dAcj) 2tt J 
2v 2 . 



-cos2A<W d4, kPcl {cf, k ) dAfaf a (Afa,4> k )cos2(Afa + j> k ). (35) 



d<j) k p c i(4>k)N a (cj) k ) 

2vr 

Realizing that elliptic flow parameter of particles from clusters is given by Eq. (TTB"]) . we rewrite Eq. (f3"5)) into 

^ = ^(l + 2v2,tV2.aCO S 2AcP). (36) 

dAcj) 2ir 

From Eq. (f3Tj) and (l36|) we obtain the total background as given by 



dN bg = N^ 
dAcp 2tt 



1 + 2v 2 ,t -^ JL V2,hy + } -T- T V 2 , a COS 2A</> 

\^bg *-f N bg 



(37) 
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where 

N bg =N hy + Y,N c iN a . (38) 



The v 2 's in Eqs. (|3"Tj) . (|3"rj)) . and (|37|) include fluctuations, so they should be replaced by v/ (i>§)- The hydro-particles 
^/ (i>|) is equivalent to two-particle «2{2} because there is no non-flow effect between hydro-particle pairs; same for 
the cluster yj (v%) because there is no non-flow effect between different clusters (we consider sub-clusters to be part 
of their parent cluster). Thus Eq. ([57)) should be 

Nbg - Nbg (1 + 2v 2 , t v 2 . bg cos 2A</>) (39) 



dA<j> 2?r 
where 



v 2 ,bg = - W JL v 2 {2} hy + 2^ — — v 2 {2} a . (40) 

We note that here cluster includes single-particle (within a give pr range) cluster. Those single-particle clusters do 
not contribute to non-flow in v 2 {2} a , but they differ from single hydro-particles because they may possess different 
v 2 values. 

In principle, v 2i t should have a similar expression as Eq. (|40[) out of symmetry reason: 

Nt.hy rm . N c i_ t N t c i_ t 

v 2 ,t = - w 2 ^-v 2 {2} tM + > — — -! — v 2 {2} t , cU . (41) 

^t.tot ,i ^t.tot 

cl_t 

where Nt.hy is number of high px trigger particles from hydro-medium (i.e., background trigger particles), v 2 {2\t,hy 
is the elliptic anisotropy of those background trigger particles, N c i_t is number of clusters of type 'cl-t' containing 
at least one trigger particle, N t . c i_t is number of trigger particles per cluster, v 2 {2} t . c i_t is elliptic flow parameter 
of trigger particles from clusters, and N tttot — N t .h y + ^^N c i_ t N t . c i_t- The only difference is that trigger particles 

cljt 

are dominated by clusters (mostly jets), and those clusters are dominated by single-trigger-particle clusters; hydro- 
medium contribution to trigger particle population should be small. We note that jet-correlation functions are usually 
normalized by total number of trigger particles including those from hydro-medium background. 

If particle correlation in clusters does not vary with cluster axis relative to the reaction plane, elliptic flow of 
particles from clusters is given by Eq. (fT5|) . or 

v 2 {2} a = v 2 {2} cl (cos2A(f>} cl . (42) 

Therefore 

v 2 , bg = ^v 2 {2} hy + V^^2{2} c /(cos2A0) c; . (43) 

N b g N bg 



IV. JET-BACKGROUND v 2 IS THE REACTION PLANE v 2 



Obviously, the elliptic flow in Eq. (|4"0")) or (|4"3"1) contains not only the two-particle anisotropy relative to the reaction 
plane, but also non-flow related to angular spread of clusters. How to obtain the elliptic flow as in Eq. (|40|) or (|43l) ? 
In Section [ill we have derived Eq. (|22|) for two-particle «2{2} in the cluster approach. The quantity in the first pair 
of parentheses in r.h.s. of Eq. (|2"2"|) is elliptic flow due to correlation with respect to the reaction plane. The second 
term in the r.h.s. arises from cluster correlation. Since elliptic flow is formally defined to be relative to the reaction 
plane, the first term in r.h.s. of Eq. (|22|) may be considered as "true" elliptic flow (except flow fluctuation effect), 
^2, flow- The second term in r.h.s. of Eq. (|22[) can be considered as non-flow, v 2> non-flow! non-flow is due to correlations 
between particles from the same dijet or the same cluster. Eq. (|2"2"]) can be expressed into 

vl{2} = u| )flow + u|, non _ flow . (44) 

Comparing Eq. (22} with Eq. (140)) , we see that 

V 2 ,bg = Inflow, (45) 

i.e., the quantity in the first pair of parentheses in r.h.s. of Eq. (|22p is the v 2 parameter in Eq. J40|) that is needed in 
constructing jet-correlation background. In other words, elliptic flow parameter that should be used in jet-correlation 
flow background is the "true" two-particle elliptic flow (i.e., due to the reaction plane and including fluctuation). 
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V. MONTE CARLO CHECKS 



In this section, we verify our analytical results by Monte Carlo simulations. We generate events consisting of 
three components. One component is hydro-medium particles according to Eq. (|27p . given hydro-particles elliptic 
flow parameter v 2 ^ y and Poisson distributed number of hydro-particles with average multiplicity N^ y . The second 
component is clusters, given cluster elliptic flow parameter i> 2jC ; and fixed number of clusters N c i; each cluster is 
made of particles with Poisson multiplicity distribution with average N a and Gaussian azimuth spread around cluster 
axis with a a . The third component is trigger particles with accompanying associated particles; the trigger particle 
multiplicity is Poisson with average N t , and the elliptic flow parameter is t> 2 ,t- The associated particles are generated 
for each trigger particle by correlation function: 



N as (4> t ) 



2n<T as 



2TT(Tns(4>t) 
( 

exp 

V 



(A0) 2 



(a<1> - n + e$ t j) 



exp 



A0-7T-0(&) 



(46) 



where the near- and away-side associated particle multiplicities are Poisson with averages N ns ((pt) and N as ((j} t ), 
respectively. The Gaussian widths of the near- and away-side peaks are fixed, and the two away-side symmetric peaks 
are set equal and their separation is fixed. All parameters in the jet-correlation function of Eq. (|46p can be dependent 
on the trigger particle azimuth relative to the reaction plane, <pf 

We first verify Eq. I|22p by generating events with hydro-particles and jet-correlated particles. (We do not include 
other clusters except jet-correlations.) We use Nh y = 150, v 2 j ly = 0.05, and N t — 2, V2,t — 0.5. We use the large 
trigger particle v 2 in order to maximize the effect of non-flow. For jet-correlation function, we generate back-to-back 
dijet with N ns = 0.7, N as = 1.2, a ns — 0.4, <j as = 0.7, and 8 — (referred to as dijet model). We fix v 2 in the 
simulation, i.e., v 2 fluctuation is not included. We simulate 10 6 events and calculate = (cos2A</>y). Including 

only hydro-particles, we obtain v 2 {2}h y = 0.05005 ± 0.00009, consistent with the input. Including all simulated 
events and all particles (hydro-particles and jet-correlated particles), we obtain v 2 {2}i nc = 0.05560 ± 0.00008. Using 
triggered events (events containing at least one trigger particle) only, we obtain v 2 {2} tr i g _evt = 0.05642 ± 0.00008. 
Using triggered events but excluding one dijet at a time (i.e., using the underlying background event of the dijet) 
and repeating over all dijets in the event, we obtain v 2 {2}t g = 0.05568 ± 0.00008. We see that the background v 2 is 
as same as that obtained from inclusive events, v 2 {2}f, g = V2{2}i nc , and both are smaller than that from triggered 
events only. 

We can in fact predict the inclusive event v 2 by Eq. (l4*3"j) using the "hydro + dijet" model. The average \J (cos2A0y) 
of jet-correlated particle pairs within the same dijet is yj (cos 2A</>jj) . £( = (cos 2A<fi)j et = 0.5054 ± 0.0004. This is 
consistent with the expected value 



(cos 2A(j>) jet 



N„ 



N ns + N a , 



■ exp I 



-2^ s ) + 



Nn 



N„ 



N ai 



exp(-2cr 2 s ) cos 29 = 0.5046 



where 6 — 0. The average yl (cos 2A(pij) of pairs of particles from different dijets is 0.2516 ± 0.0008; it equals to 
v 2 {2} a jet — v 2 ^ (cos 2 Atfi) j et = 0.5 x 0.5046 = 0.2523. The average yj (cos 2Acj>ij) for cross-talk pairs of background 
particle and jet-correlated particle is 0.10064±0. 00004; and it equals to the expected value yj v 2 {2}h y v 2t t (cos 2 A0) Jet = 
y/0.05 x 0.5 x 0.5046 = 0.1123. The inclusive event two-particle elliptic flow parameter is 



M2} 



' 150 
153.8 



x 0.05 



3.8 
153.6 



0.2523 | 



2 x 1.9 2 
153.8 2 



(0.5054 2 - 0.2523 2 ) = 0.05553; 



this is indeed consistent with v 2 {2}i nc or v 2 {2}i, g obtained from simulation. 

We now verify Eq. (j40[) or (|43]) as the correct v 2 to be used for jet-correlation background subtraction. We generate 
Poisson distributed hydro-particles with average multiplicity Ni iy = 150 and fixed elliptic flow parameter v 2 j ly = 0.05. 
We generate Poisson distributed trigger particles with average trigger particle multiplicity N t = 2.0; we use different 
jet-correlation functions (discussed below). We also include clusters that do not have trigger particles (referred to 
as minijet clusters); the particle multiplicity per minijet cluster is Poisson distributed with average N a = 5, and the 
number of minijet clusters is fixed as N c i; we fix the cluster shape to be Gaussian with width a a = 0.5 (and average 
angular spread (cos2A0) c / = exp (— 2er 2 ) = 0.6065), and also fix the cluster elliptic flow parameter w 2jC ; = 0.20. 
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TABLE I: Monte Carlo verification of analytical results of elliptic flow parameter to be used in jet-correlation background. 
Hydro-particle multiplicity, trigger particle multiplicity, jet-correlated near- and away-side multiplicities, and particle multi- 
plicity per minijet cluster are all generated with Poisson distributions, with averages N ny , Nt, N n3 , N aa , and N a , respectively, 
while number of minijet clusters is fixed as N c i. The jet-correlation function is given by Eq. (|46p . with near- and away-side 
Gaussian width fixed to be a na = 0.4 and <j aa = 0.7, respectively. The minijet cluster Gaussian width is fixed to a a = 0.5. The 
elliptic flow parameters for hydro-particles, trigger particles, and clusters are t»2,hy, V2,t, and V2 >c i, respectively, and are fixed 
over all events without fluctuation. We use Nh y = 150, Nt — 2, iV a = 5, V2,h y = 0.05, and V2 tC i = 0.20. Four cases are studied; 
the parameters for the cases are, respectively, (i) N c i = 0, V2,t = 0.5, C — 0, N n3 = 0.7, N as = 1.2, a na = 0.4, a as = 0.7, and 
9 = 0; (ii) N cl = 10, v 2 .t = 0.5, C = 0, N ns = 0.7, N as = 1.2, a ns = 0.4, a as = 0.7, and 9 = 0; (iii) N cl = 10, v 2 ,t = 0.5, C = 0, 
N ns = 0.7, N as = 1.2, a na = 0.4, a aa = 0.7, and = 1; and (iv) N c , = 10, v 2 ,t = 0.1, C* = 0, N na = 0.7, N aa = 1.2, a na = 0.4, 
a aa = 0.7, and 9 = 1. 



Case 

(i) hydro + dijet 

(ii) hydro + minijet + dijet 

(iii) hydro + minijet + near-side + 
away-side double-peak 

(iv) hydro + minijet + near-side + 
reaction-plane dependent away-side 
double-peak + clusters 



V2{2} i>2,at Calculated V2,b a 

0.05557(8) 0.05505(8) 0.05500 

0.08465(6) 0.07115(8) 0.07126 

0.08172(6) 0.06815(8) 0.06813 

0.08279 0.06883(35) 0.06910 



We simulate 10 6 events and form raw correlation functions normalized by the number of trigger particles. In order 
to extract the real background v 2 from the simulations, we subtract the input jet-correlation function. If the jet- 
correlation function varies with the trigger particle angle relative to the reaction plane, the trigger multiplicity weighted 
average jet-correlation function is subtracted. We fit the resultant background function to B (1 4- 2w 2jt W2,fit cos 2A0) 
where B and V2fix are fit parameters. We treat the input v 2 j as known; we did not include any complication into v 2> t- 
We compare the fit w 2 ,fit to the calculated one by Eq. (J40"f or (|43|) . We study several cases with different shapes for 
jet-correlation function, as well as varying values for some of the input parameters: 

(i) "hydro + dijet" model: we generate back-to-back dijets accompanying trigger particles, without other clusters. 
The calculated v 2 , bg by Eq. g3J) is v 2M = ^ x 0.05 + ^| x 0.5 x 0.5046 = 0.05500. 

(ii) "hydro + minijet + dijet" model: we include minijet clusters in addition to (i). The calculated v 2 ,bg by Eq. (|4"3l) 
is v 2 , bg = 2ufg x 0.05 + x 0.2 x 0.6065 + |g|f x 0.5 x 0.5046 = 0.07126. 

(iii) "hydro + minijet + near-side + away-side double-peak" model: we generate jet-correlated particles by correlation 
function with double-peak away-side to replicate the experimentally measured reaction-plane averaged dihadron 
correlation function [TH) We used the same Gaussian parameters for the correlation peaks as in (i) but 

= 1, thus (cos2A0) je4 = 0.1689. The calculated v 2 , bg by Eq. (gSJ) is v 2M = ^gfg x 0.05 + f^f x 0.2 x 0.6065 + 



2x1.9 
203.8 



x 0.5 x 0.1689 = 0.06813. 



(iv) "hydro + minijet + near-side + reaction-plane dependent away-side double-peak" model: we include reaction- 



plane dependent jet-correlation function similar to preliminary experimental data llj. We have to use Eq. (|40[) 
to calculate v 2 ,b g , which gives v 2M = x 0.05 + x 0.2 x 0.6065 + x 0.1423 = 0.06910. Note that, 

in this simulation of reaction-plane dependent jet-correlation signal, the number of jet-correlated particles is 
not 1.9, but rather 1.74. Also note that, due to the reaction-plane dependency of the jet-correlation signal, the 
elliptic anisotropy of jet-correlated particles cannot be factorized into the product of the trigger particle elliptic 
flow and the average angular spread of the jet-correlation signal as in Eq.[l5l but has to be calculated by Eq.[T3l 

We list our comparison in Table HI The fit v 2: at is supposed to be the real background v 2 ,bg- The fit errors are due 
to statistical fluctuations in the simulation. As can be seen, the calculated v 2 ^ g reproduces the real background v 2 ^ g 
in every case. The v 2 ,b g values differ from the hydro-background v 2 due to contributions from cross-talks between 
cluster correlation and cluster flow. Also shown in TableUare the two-particle «2{2} from all pairs in inclusive events. 
The «2{2} values differ from v 2i b g due to non-flow contributions between particles from the same dijet or the same 
cluster. 

Figure [lja) shows the raw correlation function for case (iii) and flow background using the calculated v 2 ^ g and 
normalized by ZYAl . Figure[TJb) shows the ZYAl-background subtracted jet-correlation function, using the calculated 
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v 2.bg by Eq. [40] for flow background. The background-subtracted jet-correlation is compared to the input signal. As 
shown in Fig.fjjb), the shapes of the input signal and extracted signal are the same, which is not surprising because 
the calculated v 2 ^ g is the correct value to use in flow background subtraction. The roughly constant offset is due to 
ZYAl-normalization. 




-1 1 2 3 4 5 




-1 1 2 3 4 5 




-1 1 2 3 4 5 



FIG. 1: (Color online) (a) Simulated raw correlation from the "hydro + minijet + near-side + away-side double-peak" model, 
Case (iii) in Table[I] ZYAl-normalized flow-background using the calculted V2,b g is shown as the red histogram, (b) Background- 
subtracted jet-correlation (black data points) compared to the input correlation signal (red histogram) . The background uses 
the calculated t>2,6 9 and is normalized to signal by ZYA1. The input signal shifted down by a constant is shown in the blue 
histogram, (c) As same as (b) except the subtracted background uses the decomposed V2{2D}. 



VI. HOW TO "MEASURE" JET-CORRELATION BACKGROUND (REACTION PLANE) v 2 



Two-particle angular correlation is analyzed by STAR and is decomposed into two components 12|: one is the 
azimuth quadrupole, v 2 {2D}, that is due to correlations of particles to a common source, the reaction plane; the 
other is minijet correlation that is due to angular correlation between particles from the same minijet or the same 
cluster. Properly decomposed, the azimuth quadrupole should correspond to the first term in r.h.s. of Eq. (|22[) . 

v 2 {2D} = ^v 2 {2} hy + £ ^v 2 {2} a . (47) 

This is identical to Eq. (|40p . That is, the elliptic flow parameter from a proper 2D quadrupole-minijet decomposition 
is exactly what is needed for jet-correlation background calculation. Decomposition of minijet correlation and flow, 
assuming the functional form for minijet correlation, has been carried out by STAR as a function of centrality but 
including slip? [12|. One may restrict to narrowly windows to obtain v 2 {2D} as a function of pt, however, statistics 
can quickly run out with increasing px because the 2D decomposition method requires particle pairs. 

Figure[Ifc) shows the ZYAl-background subtracted jet-correlation function, using the decomposed v 2 {2D} from the 
simulation data for flow background. The background-subtracted jet-correlation is compared to the input signal. The 
shapes of the input signal and the extracted signal are the same, which demonstrates that the decomposed v 2 {2D} is 
close to the input elliptic flow value. Again, the roughly constant offset is due to ZYAl-normalization. 

One natural question to ask is why not to decompose jet-correlation and jet-background directly from high-p^ 
triggered correlation function. One obvious reason is that jet-correlation shape is unknown a priori, thus one cannot 
simply fit triggered correlation to a given functional form. However, even when the functional form of jet-correlation 
signal is known, as is the case in our simulation, we found that the decomposed jet-correlation signal shape deviates 
significantly from the input one. This is because the jet-correlation signal is not orthogonal to flow background, but 
rather entangled, both with near- and away-side peaks, and hence one can get false minimum \ 2 m decomposing the 
two components with limited statistics. 



VII. IMPLICATIONS AND SUMMARY 



In experimental analysis, v 2 values from various methods have been used for jet-correlation background. STAR used 
the average of the reaction plane v 2 {RP} and the four-particle v 2 {4} and used the range between them (or between 
v 2 {2} and «2{4}) as systematic uncertainties [9. Il5j. The reaction plane v 2 {RP} and the two-particle v 2 {2} contain 
significant non-flow contributions, while the non-flow contributions are significantly reduced in the four particle «2{4}. 



11 



On the other hand, any v 2 fluctuation reduces «2{4}. The recently measured v 2 {2D} magnitudes from STAR are 
similar to the v 2 {4} results, suggesting that the used v 2 values for background calculation are too large, by about 
la systematic uncertainty. PHENIX used v 2 results from the reaction plane method where the reaction plane is 
determined by particles several units of pseudo-rapidity away from particles used in jet-corrclation analysis. Some 
but not all non-flow is removed; the remaining non-flow may be dominated by the long range Arj correlation (ridge) 
observed in non-peripheral heavy- ion collisions [HI, [h| [T3] - Thus the v 2 values used by PHENIX for background 
calculation are also too large. 

In summary, we have derived an analytical form for jet-correlation flow-background in a cluster approach. We argue 
that the elliptic flow v 2 parameter to be used in jet-correlation background is that from two-particle method excluding 
non-flow correlation unrelated to the reaction plane, but including cross-terms between cluster correlation and cluster 
flow. We verify our result by Monte Carlo simulation for various jet-correlation signal shapes as well as varying other 
input parameters to the simulation. We demonstrate that the v 2 parameter to use in jet-correlation flow background 
is as same as the v 2 {2D} from a proper 2D quadrupole-minijet decomposition of two-particle angular correlation. 
However, we note that 2D quadrupole-minijet decomposition requires a model for minijet correlation shape, which 
gives rise to systematic uncertainty on the extracted v 2 {2D}. 
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